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Abstract 

A scanning tunneling microscopy (STM) image of a hydrogen-adsorbed Si(OOl) surface is studied 
using first-principles electron-conduction calculation. The resultant STM image and scanning 
tunneling spectroscopy spectra are in agreement with experimental results. The contributions of 
the vr states of bare dimers to the tunnel current are markedly large, and the a states of the dimers 
rarely affect the STM images. The tunnel currents do not pass through the centers of the dimers 
but go through the edges of the dimers with local loop currents. In addition, when the tip exists 
above the hydrogen-adsorbed dimer, there are certain contributions from the vr state of the adjacing 
bare dimers to the tunnel current. This leads to the STM image in which the hydrogen-adsorbed 
dimers neighboring bare dimers look higher than those surrounded by hydrogen-adsorbed dimers. 
These results are consistent with the experimental images observed by STM. 

PACS numbers: 72.20.-1, 02.70.Bf, 68.37.Ef, 73.20.-r 



1 



There is considerable theoretical and experimental interest in the electronic and geometric 
structures of semiconductor surfaces, especially those of silicon surfaces. Thus far, scanning 
tunneling microscopy (STM) [l| has made it possible to observe the structures at surfaces 
with atomic resolution and has provided much useful information on surface configurations 
P]. It is well known that the STM image is not directly related to the atomic configuration 
but depends on the electronic structure at the surface. For example, in the case of H-, C1-, 
or NH3- adsorbed dimers of Si(OOl) surfaces, the adsorbed dimers look geometrically lower 
than bare dimers in STM images 3, 0, 0, 0] . 

On the theoretical side, first-principles calculations based on density functional theory [3] 
have been employed for identification of surface atomic structures in STM images. Analyses 
on the basis of the local density of states (LDOS) of surfaces can yield reasonable results for 
STM images under the Tersoff-Hamann (TH) 8] procedure, which means that the variation 
of the tunnel current with respect to a bias voltage is proportional to the LDOS of the surface 
at the tip position. By virtue of its simplicity, the TH approximation has been successfully 
applied to the interpretation of STM images and has provided us with reasonable results 
in many cases. However, there are some inherent uncertainties in the TH approximation. 
For example, the images generated by the TH approximation sometimes involve spatially 
localized surface states ,9] which do not contribute to electron conduction due to the discon- 
nection with the propagating Bloch states inside electrodes j,10], and there are no definite 
ways to distinguish actually dominant states to electron conduction from other localized 
states. In addition, it is difficult to examine the tunnel current path because the global 
wave functions for infinitely extended states continuing from the substrate to the tip are not 
accounted for in the framework of the TH approximation. For the accurate interpretation 
of STM images, it is mandatory to compute the tunnel current fiowing between the tip and 
surfaces. 

FIG. 1: Schematic representation of computational model. The thick solid line represents the 
supercell. Large light, small light, and dark balls are aluminum, hydrogen, and silicon atoms, 
respectively. In the top view, atoms are denoted by large and small balls according to the distance 
from the surface. 



In this Rapid Communication, the first-principles electron-conduction calculation for tun- 
nel current fiowing between the hydrogen-adsorbed Si(OOl) surface and tip of STM is imple- 
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merited. To the best of our knowledge, this is the first theoretical indication of the current 
image and current pass in STM systems using first-principles calculation. The calculated 
STM image and scanning tunneling spectroscopy (STS) spectra are consistent with exper- 
iments. We found that the surface tt states of the bare dimer largely affect the current 
image while the contributions of the a states of the dimers are quite small. In addition, 
tunnel currents flow through the edges of the dimers rather than the centers of the dimers. 
Moreover, there are certain contributions of the vr states of the neighboring bare dimer to 
the tunnel current when the tip is above the hydrogen-adsorbed dimer. 

FIG. 2: Theoretical constant-height STM image of Si(OOl) surface. In order to compare with STM 
images, positions with high tunnel current are shaded Ughter than those with low tunnel current. 
The meanings of the symbols are the same as those in Fig. ^ 



Our 

proach 



irst-pri ncip les calculation program is based on the real-space finite-difference ap- 



13|, in which boundary conditions are not constrained to be periodic and 
basis sets used to expand wave functions are not required. Therefore, one can calculate not 
only electronic structures of surfaces but also tunnel current flowing between a tip and a 
surface of STM accurately by virtue of making efficient use of the advantages in boundary 
conditions. Moreover, the wave functions in vacuum regions can be strictly determined in 
the real-space formalism, while their characteristics unphysically depend on the choice of 
basis set in the linear combination of atomic orbitals approach. Figure ^ shows the compu- 
tational model. The scattering region consists of five silicon layers with the (2x2) lateral cell 
and pyramidal tip made of aluminum atoms. We adopt the nine-point finite-difference for- 
mula for the derivative arising from the kinetic-energy operator of the Kohn-Sham equation 
in the density functional theory We take a cutoff energy of 23 Ry, which corresponds 
to a grid spacing of 0.65 a.u., and a higher cutoff energy of 210 Ry in the vicinity of nu- 
clei with the augmentation of double-grid points [l^ llS^ . The ion cores are rep resented 



by the norm- conserving pseudopotential 



14| given by Troullier and Martins llSl and the 



exchange-correlation effects are treated by the local density approximation jl6|. The su- 
percell is imposed the periodic boundary condition in the [110] and [110] directions while 
the nonperiodic boundary condition in the [001] direction. Hydrogen atoms are adsorbed 
on only one silicon dimer and the bottom side of the film remains bare. All of them are 
suspended between semi-infinite aluminum jellium electrodes. The distance between the tip 



and surface is set to be ~ 9.5 a.u. so as to save computational cost. The atomic geometry 
of the sihcon surface is determined by first-principles structural optimization using the con- 
ventional periodic supercell technique in which the bottom side of the film is terminated by 
hydrogen atoms and the tip is absence. The global wave functions for infinitely extended 
states continuing from one electrode side to the other are determined by employing the over- 
bridging boundary matching method Il3l. ll7l |. The tunnel current of the system is evaluated 



using the Landauer-Biittiker formula [ji 



2e 



Ep 



/ = -- T{E')dE', (1) 

where h is Planck's constant, e is the electron charge, Ep is the Fermi level, and T{E) is the 
transmission coefficient of incident electrons at energy E. The tunnel current is examined on 
4x4 points on the surface by displacing the aluminum tip, and the STM image is generated 
by integration of the tunnel current over the range from the Fermi level Ep to Ep — 1.5 eV. 

Figure El depicts a constant-height STM image generated by the tunnel current. The 
tunnel current at the bare dimer is larger than that at the hydrogen-adsorbed dimer. Since 
the tip must approach the sample surface in order to achieve the constant tunnel current, 
the hydrogen-adsorbed dimer looks geometrically lower than the bare dimer in STM, which 
agrees with the experimental results . It is also noted that a tunnel current of ~ 

10^^ A is consistent with the experimental result in consideration of the fact that the tunnel 
current increases by about seven times per 1 A as the tip-surface distance decreases 0,0]. 
One can recognize that the resultant image of hydrogen adsorbed dimer is not consistent with 
typical experimental image in which the hydrogen adsorbed dimer shows two protrusions 
across the dimer directions with a minimum at the center 0|. This is due to the effect of 
the adjacent bare dimer discussed later Indeed, the two protrusions look fuzzy when 
the hydrogen adsorbed dimers neighbor to bare dimers in the experimental image. 



FIG. 3: STS spectra on bare dimer (solid curve) and hydrogen-adsorbed dimer (dashed curve) in 
the computational model of Fig. ^ The zero of energy is chosen to be the Fermi level. 

In order to explore the relationship between the surface electronic structure and the 
tunnel current. Fig. El shows the calculated STS spectra when the tip is located above the 
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bare dimer and above the hydrogen-adsorbed dimer using Eq. (|T)). Considering the fact that 
the local density approximation underestimates the band gap by ~ 0.6 eV, the peaks in the 
spectra are in agreement with those in experiments ^. The vr (vr*) state of the bare dimers 
is located at about 1.5 — 0.5 (0.5 — 1.0) eV below (above) the Fermi level Therefore, we 
can conclude that the peaks in the STS spectra of the bare dimer at —0.6 and 0.6 eV are 
due to the vr and vr* states, respectively. It is noted that in the case of hydrogen-adsorbed 
dimer the peaks of tt and it* are also observed although the dangling bonds of the dimer are 
passivated. 

For further interpretation of the STM image and STS spectra, we show in Fig.|3]the charge 
distributions and the current distributions of electrons incident from the silicon surfaces to 
the tip with energies between the Fermi level Ep and Ep — 1.5 eV, and the contour plot of 
the LDOS which are obtained using the conventional three-dimensionally periodic supercell 
for the same energy range. The current distribution is derived by the equation of continuity. 
From Fig. Efa) in which electrons accumulate above and below the Si-Si bond of the bare 
dimer, one can recognize that the vr bonding state largely contributes to the tunnel current 
although there is little contribution of the a bonding state; the surface vr-states of the 
bare dimer, which are expected to be localized in the surface region, actually contribute to 
electron conduction from the sample electrode to the tip electrode in the case of the Si(OOl) 
surface. 

In addition, the local loop currents rotating around the atoms consisting the dimers are 
observed in Fig. 01^b). These peculiar characteristics of the tunnel current can be explained 
by the decomposition of the scattering wave functions into the atomic orbitals around the 
respective atoms. 



Here, x, y and z directions correspond to [110], [110] and [001] directions, respectively. By 




(2) 
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substituting Eq. (j2I) into the equation of continuity, we have 
I{r) 

eh 



r MC'^)Cl-ct>l{r)Vct>,Sr)]dE' 
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a=x,y,z 



Ef 



Je 



a,p — x,y,z 



(3) 



The first term of Eq. (jS)) corresponds to currents fiowing along the a direction and the second 
one to local loop currents rotating on the a-f3 plane. Due to the geometrical symmetry of 
atoms, the currents except that fiowing along the z direction vanish inside the substrate. On 
the other hand, the current rotating on the x-z plane emerges on the surface plane because 
the atoms on the surface plane form dimers. Consequently, the currents at the edges of the 
dimers are enhanced as seen in Fig. Efb). The wave functions of the surface can be also 
expressed in terms of the following basis set, 

0. = [0.(r) + 0p^(r)]/v^ 



= [-0,(r) + v^0,^(r) + 0,^(r)]/2 

0<x,, = i-Ur) - V2<t>p,ir) + 0p,(r)]/2. (4) 

One can notice that the vr state contributes to the current along the z direction while the 
cr state affects the local loop current around the surface atoms. This agrees with the STS 
spectra of Fig. 01 and the charge distributions of Fig. Ufa), in which the vr states actually 
contribute to the tunnel current. 

Moreover, it is observed in the lower panel of Fig. El^b) that there is a certain contribution 
from the vr state of the neighboring bare dimer to the tunnel current, which agrees with the 
STS spectra of Fig. El where the vr state affects the spectra even when the tip locates above 
the hydrogen-adsorbed dimer. In order to estimate the contribution of the neighboring 
bare dimer, we evaluate the following quantity, R = J^^ Iz{r)dS/ ^ Iz{r)dS, where the 
surface integral regarding Sbar (Sads) is implemented on the rectangle (2x1) plane just above 
the bare (hydrogen-adsorbed) dimer. The i? is ~ 0.25 around the bisection of the surface- 
dimer- and tip-apex-atom planes. When the tip-surface distance increases to 10 A, the 



FIG. 4: (Color) (a) Charge distributions and (b) current distributions, —I{r)/e, of incident elec- 
trons from the silicon surfaces to the tip. The aluminum tip is located above the bare dimer in the 
upper panel and above the hydrogen-adsorbed dimer in the lower pannel. The density distribu- 
tions of electrons incident from the bottom electrode with energies between and Ep — 1.5 eV 
are depicted, (c) Contour plots of the LDOS, which are computed using the conventional three- 
dimensionally periodic supercell, for the same energy ranges as those in (a). The channel electron 
distributions and the contour plots of the LDOS are illustrated on a logarithmic scale and each 
contour represents twice or half the density of the adjacent contour lines. The map is along the 
cross section indicated by the dotted lines of A, B, or C in Fig. ^ Gray, pink, and blue balls are 
aluminum, hydrogen, and silicon atoms, respectively. The large (small) balls indicate the atomic 
positions on (above) the cross section. 

R becomes 0.53 at the bisection plane, which implies that 35 % of the tunnel current is 
attributed to the neighboring bare dimer: The transmission from the vr states to the tip 
reduces more moderately than that from the states around the hydrogen-adsorbed dimer 
as the tip-surface distance increases because the vr state is widely expended to the vacuum 



in which 



region. In addition, our result is also consistent with the height profile of Ref. 
hydrogen- adsorbed dimers neighboring bare dimes look geometrically higher by ~ 0.25 A 
than those surrounded by hydrogen-adsorbed dimers in the [010] direction. 

Let us discuss the relationship between the charge distribution of the incident electrons 
and the LDOS. In the case of the bare dimer, there is no significant difference in charac- 
teristics of the counter plot of LDOS in Fig. |31^c) from that of the charge distribution of 
incident electrons in Fig. Ufa). Although a certain influence of the adsorbed hydrogen atoms 
on the LDOS is observed around the Si-H bonds in Fig. IHc), the contributions of the states 
affected by the hydrogen atoms to the distribution of incident electrons are small compared 
with the Si-Si bonding state in Fig. 131(a); the densities around the hydrogen atoms and Si-H 
bonds are one forth of that around the Si-Si bond in Fig. El(a) while that is half in Fig. |31^c). 
We assured that this situation does not change in the case of large tip-substrate distance. 
This result implies that by being adsorbed hydrogen atoms the orbitals comprised by the 
s and Pz orbitals of the surface atoms couple weakly with the waves continuing from the 
substrate compared with those composed by the Px orbital. 
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In summary, we have studied an STM image of hydrogen-adsorbed Si(OOl) surfaces us- 
ing first-principles electron-conduction calculation. The tunnel current flowing between a 
hydrogen- adsorbed dimer and the STM tip is smaller than that flowing between a bare 
dimcr and the tip. The a states of the dimer hardly affect the STM image while the con- 
tributions of the 7r bonding states to the tunnel current between the bare dimer and the 
tip are large. In addition, the currents passing through the edges of the dimers are due to 
the contributions of the a states of the dimers. Since the tunnel current with the tip being 
above the hydrogen- adsorbed dimer is attributed to the tt states of the adjacent bare dimers, 
hydrogen- adsorbed dimers neighboring bare dimers appear geometrically higher than those 
surrounded by hydrogen- adsorbed dimers in the [010] direction. These results are consistent 
with experiments. 

This research was partially supported by a Grant- in- Aid for the 21st Century COE "Cen- 
ter for Atomistic Fabrication Technology" , by a Grant-in- Aid for Scientific Research in Pri- 
ority Areas "Development of New Quantum Simulators and Quantum Design" (Grant No. 
17064012) and also by a Grant-in-Aid for Young Scientists (B) (Grant No. 17710074) from 
the Ministry of Education, Culture, Sports, Science and Technology. The numerical calcu- 
lation was carried out by the computer facilities at the Institute for Solid State Physics at 
the University of Tokyo and the Information Synergy Center at Tohoku University. 



[1] G. Binnig et al, Phys. Rev. Lett. 49, 57 (1982). 

[2] D. Drakova, Rep. Prog. Phys. 64, 205 (2001); Werner A. Hofer et al, Rev. Mod. Phys. 75, 

1287 (2003), and references therein. 
[3] R.J. Hamers et a/., Phys. Rev. Lett. 59, 2071 (1987); P. Avouris et al, J. Vac. Sci. Technol. B 

5, 1387 (1987); R.J. Hamers et al, J. Vac. Sci. Technol. A 6, 508 (1988); J.J. Boland, Phys. 

Rev. Lett. 67, 1539 (1991); J. Vac. Sci. Technol. A 10, 2458 (1992); Science 262, 1703 (1993); 

Adv. Phys. 42, 129 (1993); M.J. Bronikowski and R.J. Hamers, J. Vac. Sci. Technol. A 13, 

777 (1995); I. Lyubinetsky et al, Phys. Rev. B 58, 7950 (1998); M.Z. Hossain et al, Phys. 

Rev. B 68, 235322 (2003). 
[4] Y. Wang et al, J. Vac. Sci. Technol. A 12, 2051 (1994). 
[5] D. Chen and J.J. Boland, Phys. Rev. B 65, 165336 (2002). 



8 



[6] Z. Hu et al, Phys. Rev. B 68, 155418 (2003). 

[7] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

[8] J. TersofFand D.R. Hamann, Phys. Rev. B 31, 805 (1985). 

[9] We regard here the states spatially localized around the surface as surface states even if they 

do not exit in energy band gap in the whole Brillouin zone. 
[10] C. Noguera, J. Microsc. 152, 3 (1988); Phys. Rev. B 42, 1629 (1990); K. Makoshi, Jpn. J. 

Appl. Phys. 33, 3657 (1994); K. Kobayashi, Phys. Rev. B 65, 035419 (2002). 
[11] J.R. Chelikowsky et al, Phys. Rev. Lett. 72, 1240 (1994); Phys. Rev. B 50, 11355 (1994). 
[12] T. Ono and K. Hirosc, Phys. Rev. Lett. 82, 5016 (1999); K. Hirose and T. Ono, Phys. Rev. 

B 64, 085105 (2001); T. Ono and K. Hirose, ibid. 72, 085105 (2005); ibid. 72, 085115 (2005). 
[13] K. Hirose, T. Ono, Y. Fujimoto, and S. Tsukamoto, First- Principles Calculations in Real-Space 

Formalism, Electronic Configurations and Transport Properties of Nanostructures (Imperial 

College Press, London, 2005). 
[14] We used the norm-conserving pseudopotentials NCPS97 constructed by K. Kobayashi. See K. 

Kobayashi, Comput. Mater. Sci. 14, 72 (1999). 
[15] N. Troullier and J.L. Martins, Phys. Rev. B 43, 1993 (1991). 
[16] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

[17] Y. Fujimoto and K. Hirosc, Phys. Rev. B 67, 195315 (2003); T. Ono and K. Hirose, Phys. 

Rev. B 70, 033403 (2004). 
[18] M. Biittiker et al, Phys. Rev. B 31, 6207 (1985). 

[19] The usage of small supercell results in the discrepancy between the theoretical and experi- 
mental. When a plain TH approximation is employed using the supercell with the same lateral 
length, two clear protrusions can not be observed. If we adopt larger supercell or the surfaces 
which are fully hydrogenated, we can obtain the images which are consistent with experimen- 
tal ones. Since our interest here is the tunnel current from the surface in which both bare and 
hydrogenated dimers exist, we employ the model where one dimer is bare and the other is 
hydrogenated. 



9 



This figure "figl.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0603605vl 



This figure "fig2.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0603605vl 



This figure "figS.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0603605vl 



This figure "fig4.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0603605vl 



